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Abstract The opening profile of the cracks produced in the Double Cleavage Drilled 
Compression (DCDC) specimens for brittle materials is investigated. The study is 
achieved by combining Finite Element simulations of a DCDC linear elastic medium 
with experimental measurements by crack opening interferometry on pure silica glass 
samples. We show that the shape of the crack can be described by a simple expression 
as a function of the geometrical parameters of the sample and the external loading 
conditions. This result can be used to measure accurately in real time relevant quanti- 
ties during DCDC experiments, such as the crack length or the stress intensity factor 
applied to the specimen. 
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1 Introduction 

The Double Cleavage Drilled Compression (DCDC) sample refers to a parallelepipedic 
column with a circular hole drilled through its center that is subjected to axial com- 
pression. This geometry provides many advantages in studying the fracture of brittle 
materials. Janssen originally introduced such a fracture test for measuring the fracture 
toughness of glass (Janssen 1974). Under a uniform axial compression, the Poisson 
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effect produces a tensile stress concentrated around the central hole. This induces the 
initiation of two symmetric mode I cracks generated at each crown of the hole and 
propagating along the mid-plane of the sample as the axial compression is increased 
(see Figs. [T] and [7]). The high stability of the crack propagation has made of the DCDC 
a rather popular fracture test geometry used in various contexts including studies on 
sub-critical crack propagation in oxide glasses (Celarie et al. 2003; Bonamy et al. 2006; 
Grimaldi et al. 2008; Fett et al. 2008), on crack healing in polymers (Plaisted et al. 
2006) or on failure behavior under mixed mode loading for linear elastic materials 
(Larnder et al. 2001; Fett el al. 2005). 

The present study presents a coupled theoretical and experimental investigation of 
the crack opening profile in the DCDC specimen. A 2D Finite Element (FE) simula- 
tion based on linear elasticity is compared with a direct measurement by crack opening 
interferometry on a silica glass specimen. The results are expressed in a simple form as 
a function of the geometrical parameters of the specimen following a Williams expan- 
sion series. This experimental technique, combined with the results of the simulations, 
allows the direct measurement of the stress intensity factor driving the crack propaga- 
tion in the DCDC specimen. In addition, it provides a solid base for the interpretation 
of recent experiments on the nanoscale capillary condensation of water inside sharp 
cracks in glass (Grimaldi et al. 2008). We should stress that this linear elastic mod- 
elling only applies for very brittle materials like glass (for which the DCDC sample 
is mostly adapted). When dealing with more compliant materials such as PMMA the 
second order terms due to rotation become very relevant and should be considered by 
a nonlinear modelling (Plaisted et al. 2006). 



2 Finite Element simulation of DCDC specimens 

2.1 Loading configuration and finite element mesh 

The finite element simulation used to estimate the crack opening profile of the DCDC 
specimens is based on a 2D plane strain analysis implemented by the finite element code 
CASTSlvQ. The mesh is constituted of linear elastic isotropic four nodes elements. The 
plane strain condition is chosen to model the behavior of the sample in its mid-plane, 
where the experimental measurements are made. The details of the sample loading 
configuration and the mesh geometry are shown in Fig. [1] The sample consists of a 
prism of dimensions 2w x2tx2L with a cylindrical cross hole of radius R drilled through 
the specimen (thickness 2t). The sample is loaded with a compressive stress a applied 
to the two opposite faces. The Poisson's ratio and the Young's modulus were set to 
v = 0.17 and E = 72 GPa respectively to represent the properties of silica glass used 
for the experiments. However, the influence of E and u on the simulated displacement 
fields and the stress intensity factor Kj are explicitly given in the following. There 
are effects of v on the computed value of Kj that are not captured by the analytical 
description proposed, but they remain inferior to 1% in the range 0.1 < < 0.4. During 
the test, two symmetric cracks propagate from the central hole in opposite directions 
along the midplane of the sample. The difference in the lengths of the two cracks is 
less than 1% and we will call a the average crack length as measured from the side 



^ CAST3M software i s developed by CEA-Saclay, France. Reference web page: 
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of the hole. Thanks to symmetry of the geometry, one can hmit the simulation to one 
quarter of the specimen (composed of about 3500 nodes) , imposing specific boundary 
conditions to the system (cf. Fig. [T] bottom), i.e. Ux = along the symmetry plan 
normal to ex and Uy — along the one normal to ey on the uncracked ligament. No 
displacement constraint is imposed on the stress application faces, meaning that friction 
at the platen/sample surfaces is neglected. The presence of frictional confinement might 
partially suppress the Kj rise for cracks approaching the end of the sample (cf. Fig.[2]a), 
but the rapid acceleration of the crack at the end of the tests is indeed experimentally 
observed. In any case, the friction effect only afi'ects very long cracks that are out of 
the domain retained here. The gap between the millimeter length scale of the external 
mechanical loading and the nanometric length scale of the crack opening near the 
tip requires an adapted mesh for the calculation of stress and displacement fields. To 
analyze such a multi-scale problem without exceeding reasonable computational times, 
the elements of the mesh are chosen so that their size decreases exponentially while 
approaching the crack tip, down to a minimum element size of 0.01 nm. This allows 
reproducing correctly the asymptotic square root crack opening displacement profile 
predicted by linear elasticity down to a nanometer distance from the crack tip. 




front 



Fig. 1 Sketch of the DCDC specimen and details of the numerical simulation: The mesh used 
in the finite element calculations is only one quarter of the whole sample represented in grey 
on the scheme of the geometry. The element size decreases exponentially when approaching 
the crack tip, so that deformations of the crack profile down to one nanometer from the crack 
tip can be computed 



To investigate the influence of the geometrical parameters on the crack opening 
profile, a first series of meshes was designed with seven values for the hole radius R 
ranging from 0.4 mm {w/R = 5) to 0.8 mm {w/R = 2.5), keeping the width and length 
constant with w = 2 mm and L = 20 mm. The effect of the specimen width on the 
opening profile is studied using another series of seven meshes, varying w from 1 mm 
to 4 mm while keeping radius and length constant with R = 0.8 mm and L = 20 mm. 
For each new geometry, the crack length a was changed between 2.5 and 18.0 mm with 
0.1 mm steps. 

2.2 Evaluation of the stress intensity factor 

Once the stress, strain and displacement fields are calculated for the mesh correspond- 
ing to a given crack length a, the elastic displacement Uy{X) = Uy{x — —X,y — 0) 
{x < 0) along the crack lip is measured as a function of the distance X from the crack 
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tip. Let us note that with these notations, the full opening between the two crack lips 
is given by 2uy. The stress intensity factor Kj applied to the specimen is calculated 
using two different methods. At first, the crack opening profile is fitted using the Irwin 
equation (Ewalds 1985) 



Ki 8X 

UyiX) = —J — 



(1) 



with E' = E/{1 — v^) due to plane strain condition. The regression is done in the 
very vicinity of the crack tip, typically in the range 1 nm < X < 1 fim where the 
crack profile follows a square root behavior as shown in the following. On the other 
hand, the J-integral method (Rice 1968) is computed, providing the mechanical energy 
release rate Gj related to the stress intensity factor by the expression Gj = Kf/E' 
(Ewalds 1985). The obtained values normalized by a^/nR axe represented in Fig. [2]a 
as a function of the normalized crack length a/R for different values of the ratio w/R. 
Both methods are found in agreement within 0.5 %. 

Following He et al. (1995) and Fett et al. (2005), the dimensionless stress intensity 
factor can be expressed as 
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The previous expression was used in He et al. (1995) to fit the simulations on the domain 

2 < w/R < 4 obtaining the parameter set dg = 0, di = 1, d2 — 0.235, and d^ = —0.259. 
In Fett et al. (2005), the results were fitted on the same domain 2 < w/R < 4, but 
taking also in consideration the effect of an offset 6 of the hole with respect to the center 
of the sample {b — here) and obtaining do = —0.3703, di = 1.1163, d2 = 0.2160, 
and d3 = —0.1575. These two equations are plotted respectively in green and red in 
Fig. [2]b together with the results of our simulation. While the agreement is correct for 
some specific values of w/R (close to 4 for the expression of He et al. (1995), closer to 

3 for Fett et al. (2005)), some differences are observed, especially when reaching the 
extreme values of w/R. 

This difference is attributed to the fact that Eq. ([2} does not provide an accurate 
fit on the whole parameter range. However, an excellent fit can be obtained for a larger 
parameter range 2.5 < w/R < 5 by allowing for a second order dependence in w/R 
using the following expression 
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which leads in our dataset to the parameters cq — 0.3156, ci = 0.7350, C2 = 0.0346, 
C3 — —0.4093, C4 — 0.3794, and C5 — —0.0257. The agreement is excellent as shown by 
the comparison with the results of the simulations (blue curves) presented in Fig. [2]b. 

The deviations from the fit for the largest values of a/R are caused by the interac- 
tion of the crack with the specimen boundary and are beyond the domain of validity 
of the equation, which is given hy w < a < L~ 2w. This regime with rapid increase of 
the stress intensity factor with crack extension corresponds to the sudden and instable 
crack propagation observed in the experiments prior to complete failure of the DCDC 
specimens in two symmetric pieces. 
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Fig. 2 (a) Comparison of the values of Kj obtained by FE simulation from the crack opening 
profile (blue) using Eq. JT} and the J-integral method (red); (b) Comparison of the values of 
Kj obtained in our simulations (blue) with the formula proposed here in Eq. |3] (black) and 
the ones proposed in He et al. (1995) (green) and Fett et al. (2005) (red) 



2.3 Williams expansion series of the crack opening profile 

The crack opening profile uy{X) is first analyzed in the same geometry as the samples 
used in the experimental part of this study. In particular, we choose R — 0.531 mm, 
w — 2 mm and L — 20 mm, with an external applied stress a = 100 MPa. The 
crack opening profile obtained for a crack length a — 5 mm is shown in Fig. [S] From 
the blowup (c), we can observe that the 1 % range of validity of the Irwin equation 
(blue dashed line) is limited to the first 20 fim near the crack tip. In order to provide 
an expression with a larger range of validity, Uy(X) can be expressed in terms of a 
Williams expansion (Williams 1957; Maugis 1999) 

(»)/v-^ \^l/2 , v-3/2 , v-5/2 , , vn/2 ^.x 

Uy{X)^piX' +P3X ' +P5X ' +...+pnX' (4) 
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Fig. 3 Comparison between the crack opening profile obtained by FE simulation and the 
Williams series of order 1, 3, and 5 at different scales: (a) at the fim scale; (b) at the millimeter 
scale; (c) on the whole length of the crack. The vertical lines at 15 fim, 300 ^m and 1.7 mm 



represent the 1 % limit of validity of u. 
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where only the odd terms pi, p3, ps,... pn are present due to the symmetry conditions of 
the pure mode I loading (the displacement function Uy{x,y) must be antisymmetric in 
relation to the transformation y — > —y). The first term of this development corresponds 
to the Irwin expression given in Eq. ((TJ. As a result, the first coefficient pi is directly 
associated with the singular stress behavior and can be expressed in terms of the stress 
intensity factor 

E' 



Pi 



(5) 



As shown in the previous section, Kj (resp. pi) can be expressed as a function of the 
external stress and the geometrical parameters using the Eq. ([3]). 

The other coefficients entering in the development of the crack profile are obtained 
by an iterative fit procedure: At first, the development to the first order is sub- 
tracted to the crack profile Uy and is fitted by a function varying as X^^"^, providing 
the value of the coefficient P3. The procedure is then iterated fitting Uy — u^y^ — u^^ 
by a function evolving as X^^^, providing the value of the coefficient P5. At this fifth 
order, the agreement with the results of the simulation is excellent up to X ~ 1.7 mm 
as shown in Fig. [3]c. As a result, we limit our analysis to this order in the following of 
this study. The limitations of the developments of lower orders are presented in Fig. [3]a 
and[3]b. In particular, the Irwin equation - development to the first order - is found 
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to be valid for X < 15 fim. The range of validity of these Williams developments is 
not universal and does depend on the sample dimensions as discussed in section ID 




Fig. 4 Crack opening profiles obtained for different crack lengths. Once normalized by pi, 
proportional to the stress intensity factor (see Eq. l|5j), all the profiles follow the same master 
curve 



The effect of the crack length on the crack opening profile is now investigated. Fig. 
U] represents the profile obtained for five values of a ranging from 5 mm to 15 mm. As 
shown in the inset of Fig. [J] all profiles have the same shape once normalized by the 

first coefficient pi = ^^\J~^ of their Williams development. In other words, the crack 
opening profiles can be written as 

'^^(^^ = (^'^' + + "^^'^') 

where the coefficients 03 = ^ and cib ~ ^ are independant of the crack length. 

Very interestingly, these coefficients are also independent of the radius of the hole 
in the DCDC specimens. This is shown in Fig. |5] where the variations of 0:3 and 0:5 
with the crack length are plotted for seven specimen geometries with different radius 
R, keeping constant the specimen width w = 2 mm. All the curves present an extended 
plateau regime as far as the crack tip is not too close to the specimen sides, as expected 
from the previous observations. Moreover, all the plateau regimes correspond to the 
same value of 0:3 and as, irrespective of the radius of the hole. In other words, hole 
radius and crack length have no infiuence on the shape of the crack lips in a broad 
range of these geometrical parameters. 

However, a dimensional analysis suggests that both coefficients must vary with a 
length of the problem. More precisely, 01^ and 05 are expected to scale as the inverse 
of a length and the inverse of the square of a length, respectively. The specimen width 
being the only remaining relevant length scale of the problem, it is natural to observe 
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Fig. 5 Normalized coefficients 03 and 0^5 involved in the Williams development of the crack 
opening profile (see Eq. If the crack tip is sufficiently far from the specimen sides, these 

coefficients, and therefore the shape of the opening profile, are independent of both the crack 
length a and the radius R of the hole of the DCDC specimens 



on Fig. [6] the following scalings 

1.319 , 0.515 
as = and Q5 = ^ (7) 

w 

where the dimensionless constants are provided by linear regressions represented by 
straight lines. 

As a result, combining Eq. (|6]) with Eq. ((T}, one gets a rather simple description of 
the crack opening profile 



.,(X)^|fy|vx(l + 1.319g)+0.515g)^). (8) 
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Fig. 6 Variations of the normalized coefficients 03 and 05 involved in the Williams develop- 
ment of the crack opening profile (see Eq. l[6]l) with the specimen width w. The straight lines 
corresponds to variations as aa ~ 1/to and 05 ~ as given in Eq. Q 



This expression remains valid with a precision better than 1 % as far as the effects of 
the finite size of the specimen are not relevant, i.e. for 2.5 w — R < a < L— 2.3 w — R, and 
relatively close to the crack tip, i.e. for X < 0.85 w, because we limited the Williams 
development to the fifth order. 

Using the Eq. (|3]) providing the expression of the stress intensity factor as a function 
of the geometrical parameters, one obtains here a closed form of the opening profile 
as a function of simply measurable quantities of the problems. This can be very use- 
ful during DCDC experiments. Indeed, combining this analytical expression with the 
experimental measurement of the crack opening profile using e.g. the interferometry 
technique described in the next section can provide an easy way to assess accurately 
the stress intensity factor that drives the crack growth in DCDC experiments. 

But before discussing in more details these applications, we compare now the main 
results of this analysis with experimental measurements. We focus on one single DCDC 
geometry and compare the opening profile measured experimentally with the result of 
the finite element simulations. 



3 Experimental measurement by crack opening interferometry 

The experimental measurement of the crack opening profile is performed by using a 
precision loading apparatus (based on a Microtest load cell produced by Deben, Wool- 
pit, UK) providing highly stable failure conditions necessary for in-situ observations 
of slow crack propagations in glasses by Atomic Force Microscopy (AFM) (Celarie 
2004). DCDC samples of pure silica glass (Suprasil 311, Heraeus, Germany) are ma- 
chined into 4 X 4 X 40 mm'^ bars {w = 2 mm and L = 20 mm) and polished with 
Ce02 to a RMS roughness of 0.25 nm for an area of 10 x 10 fim^. A hole of radius 
R = (0.531 ±0.010) mm was drilled at their center to trigger the start of two symmetric 
fractures (see Fig. [T}. 

The measurement of the crack opening profile is performed by means of monochro- 
matic light reflection interferometry (Sommer 1970; Liechti 1993) according to the 
sketch of Fig. [S] When observing the crack tip by AFM at the free surface of the sam- 
ple, the crack front is vertical as in Fig.[7]b. In order to use the AFM setup to measure 
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Fig. 7 Experimental setup: (a) Sketch of the DCDC sample arrangement in the crack opening 
interferometry technique (horizontal crack plane); (b) picture of the experiment in the AFM 
in-situ crack tip imaging mode (vertical crack plane) 



the crack opening displacement, the DCDC sample had to be placed with the frac- 
ture front lying horizontally. The white light source (illuminating vertically the sample 
after passing through the optical axis on the CCD camera and being reflected by an 
inclined mirror in the AFM head) was substituted by a green laser source (wavelength 
A — (532 ± 1) nm) in order to provide an intense monochromatic light (Fig. [7] a). 




Fig. 8 Sketch of the crack opening interferometry technique. NB: exaggerated vertical scale 
for the crack opening and lateral scale for the reflected beams 



When the monochromatic beam reaches the fracture plane in an orthogonal direc- 
tion, it is split into two beams that are reflected by the two fracture surfaces and then 
collected back by the CCD camera. The path difference between the two beams is twice 
the distance between the two lips of the crack, i.e. A = 2d = 4uy{X) according to the 
notations used previously. The phase difference 5 is therefore given by 

S^n + ^4uy{X) (9) 

where n is the refractive index of the fluid flUing the gap in the fracture and the term 
TT is a constant additive phase shift caused by the fact that both beams are normally 
reflected on symmetrically opposite interfaces (respectively glass/fluid and fluid/glass) 
(Born and Wolf 1999). 

The intensity I of the reflected light is given by: 
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I^h+h + 2{hh)^'^ cos(5) ~ 4Ji cos^ (0 (10) 

where the approximation on the right stems from the fact that the two reflected beams 
have similar intensity. Since Uy{X) is a monotonically increasing function of X, the 
reflected intensity will develop into a series of fringes parallel to the crack front, each 
new one corresponding to an increase of the crack opening d = 2uy{X) by a quantity 
A/2. The term n in Eq. @ implies the presence of an intensity minimum at X = 0. 
The first intensity maximum is obtained for 2uy{Xi) = A/4 and the position X^, of 
the fc*'' order fringe is given by the equation 

2«y(^fe) = (i + fc) ^. (11) 

By measuring the position and order of the fringes along the whole crack, we obtain 
a series of discrete data points on the full crack opening profile as shown in Fig. 15] 

In order to provide a better resolution, the fringe pattern along the 5 mm crack 
length is imaged with a series of 1 mm overlapping images. To insure that we preserve 
the right order and position of the fringes, the displacements of the load cell are mea- 
sured by a heterodyne interferometer (Zygo) with 10 nm resolution. The uncertainty 
in the position of each fringe (~ 50 fim) was estimated depending on the quality of the 
local signal and fringe shape (the error bars are plotted in Fig. [O]). A special attention 
must be paid to the first fringe, since the opening gradient is maximal in that region. 
However, in our operating conditions, the position of the first fringe is located to some 
hundreds of micrometers from the crack tip, making its measurement safe. The even- 
tual optical bias induced by the stress gradient on the apparent crack tip position can 
be shown to be less than 0.1 ^m in glass samples (Kysar 1998). 




X (mm) 

Fig. 9 Comparison between the opening crack profile obtained by the four experiments (in 
color, with errorbars) at different loads (renormalized) and the results of the FE simulation 
for the same set of parameters 

In order to verify the robustness of the method and the linearity of the loading 
configuration and sample, the measurement was repeated for four different loads (a — 
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50, 62.5, 75 and 87.5 MPa) and the opening profile was normalized to the reference load 
a = 100 MPa also used in the FE simulations according to linear elasticity. The four 
experimental series were realized in a time span of one hour in order to guarantee stable 
conditions (T = (21 ± 1)° C, RH — (35 ± 5) %) and minimal crack propagation. Since 

1 /2 

the highest stress used in this experiment corresponds to Ki — 0.36 MPam ' , the 
maximum crack propagation velocity is of the order of 10"^" ms^^ (as determined by 
AFM in-situ crack propagation experiments effectuated on the same sample as shown 
in Fig. [7]b, cf. Celarie 2004). The maximum crack advance in one hour can thus be 
estimated to 0.4 /im, which is below the resolution of our optical microscope. The crack 
length can thus be considered as constant with the value a — (5.02 ± 0.01) mm. 

The results of Fig. |9] show an excellent agreement between the four measurements 
within the experimental errors. The agreement is also excellent when comparing the 
experimental measurements to the results of the FE simulation for the same set of 
parameters (represented by a black line on Fig. [9]). We point out that the agreement 
with the FE simulations is also excellent in the diverging part near the central hole 
(corresponding to distances X ~ a — 5 mm from the crack tip) . We should remind that 
the FE simulation is made in 2D under plane strain condition in order to represent 
the deformation field in the bulk of the DCDC specimen where the crack opening 
interferometry is applied. When we refer to the crack length a in a 3D specimen, this 
should thus be measured in the center of the specimen. 

Finally, due to the robustness of the crack opening profile shown in the previous 
section, the agreement between the FE simulation and the experimental measurements 
for a given sample geometry can be considered as representative for a broad set of 
geometrical parameters. 

4 Discussion 

The expression ((8} for the crack opening profile in the DCDC specimen shown in 
the theoretical part of this paper and validated in the experimental section is quite 
interesting. The first important remark is that the fifth order Williams expansion is 
invariant with respect to variations in the crack length and the hole radius. Once the 
stress intensity factor is determined by Eq. this acts simply as a multiplicative 
factor to a constant crack opening profile. Of course, this is only valid for sufficiently 
long crack lengths (a > 2.5 w — R) and crack tips far enough from the specimen side 
(a < L — 2.3 w — R). In addition, this expression gives an accurate description of 
the opening profile in the domain of validity of the fifth order Williams expansion 
X < 0.85 w. 

As a consequence, the ranges of validity of the first and third order developments are 
also invariant for samples of same width w, and are expected to scale with w. The range 
of 1 % validity of the first order term (Irwin equation ((l|) is given by X < 0.0076 w, 
while the range of validity of the third order development is X < 0.15 ui. An estimate 
of these upper bounds can be easily found comparing the terms of the expansion of the 
crack opening profile, as given in Eq. ([SJ. 

We should point out that the use of typical optical systems available in AFM in- 
struments does not allow the interferometric measurement of the close neighborhood 
of the crack tip. The custom coupling of this system with a high quality optical and 
illumination system would allow using higher magnifications and extracting valuable 
information from the interpolation of the fringe intensity profiles (cf. Swadener and 
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Liechti 1998). However, the present results allow us to measure the stress intensity 
factor by a simple and direct experimental measurement of the crack opening displace- 
ment over the first millimeters from the crack tip. In fact, fitting the experimental 
profile by the Eq. (|8]) using two free parameters (position of the crack tip and multi- 
plicative constant of the profile), this technique enables an accurate measurement of 
the stress intensity factor and the crack length. This single measurement might be an 
interesting alternative to the load measurement combined with the optical estimation 
of the crack length necessary to assess the stress intensity factor, as classically done for 
DCDC fracture tests. This feature is particularly interesting for in-situ investigations of 
the crack tip since it allows reducing the weight and dimensions of the loading fixture. 

We remind that, as discussed by Plaisted et al. (2006), linear modelling can not 
provide a complete description of the DCDC deformation field, since the second order 
terms due to rotation can become very relevant, especially for long cracks. However, 
these effects are relatively limited in very stiff and brittle materials and do not seem to 
affect the present analysis as confirmed by the excellent agreement between the linear 
analysis and the experiments realized at four different loads. Moreover, subcritical 
propagation K(v) curves performed repeatedly on different crack length regions of the 
same DCDC glass specimen are found to be in satisfactory agreement between them 
(Celarie 2004), thus confirming the validity of the present treatment for glasses. A 
generalization of our approach to softer materials by accounting for nonlinear terms 
would be an interesting future development. 

Finally, it is interesting to discuss the limitations of our approach at much smaller 
scales. Due to the elevated degree of homogeneity and brittleness of silica glass, the 
lower limit of validity of the Irwin relation for a sharp crack is expected to be of 
the order of a few nanometers from the crack tip. In addition, the FE simulation is 
strictly a good representation of a measurement of the crack opening in vacuum. In 
ambient condition, capillary condensation is expected to develop into the crack cavity 
(Wondraczek et al. 2006; Grimaldi et al. 2008; Ciccotti et al. 2008). On the one hand, 
this can bias the apparent position of the crack front due to the weak reflectivity of 
the air/glass interface. On the other hand, this can alter locally the crack opening 
profile due to the attractive capillary forces. For the experiments presented here, the 
largest length of the capillary condensed phase can be estimated to few micrometers. 
As a consequence, this corresponds to a distortion of the crack lip on the first few 
micrometers close to the tip, which can unfortunately not be measured by the present 
setup, but which could be revealed by the coupling it with a higher quality optical and 
illumination system. The effect of capillary condensed phase represents an interesting 
challenge for future experimental and theoretical investigations. 

5 Conclusion 

This work reports a combined theoretical and experimental study of the crack opening 
profile in DCDC specimens of linear elastic brittle materials. It has been shown how a 
fifth order Williams expansion series can be used to describe very accurately the crack 
opening profile over a conveniently large domain near the crack tip. Finite element sim- 
ulations have allowed on one side to obtain an accurate expression of the variations of 
the stress intensity factor with the crack length for a broad range of DCDC geometries, 
and on the other side, to express the crack opening profile in a simple form as a function 
of the geometrical parameters of the specimen. The ranges of validity of these analyti- 
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cal expressions and more specifically of the Williams development at various orders are 
calculated and interpreted. Moreover, the simulated profiles were successfully tested 
against an accurate experimental measurement by monochromatic-light crack opening 
interferometry. This single optical measurement might be an interesting alternative to 
the load measurement combined with the optical estimation of the crack length both 
necessary to assess the stress intensity factor, as classically done for DCDC fracture 
tests. 
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